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Our local environment at r < 10 Mpc expands linearly and smoothly, as if 
ruled by a uniform matter distribution, while observations show the very clumpy 
local universe. This is a long standing enigma in cosmology. We argue that the 
recently discovered vacuum or quintessence (dark energy (DE) component with the 
equation of state pq — wpQC 2 , w G [—1,0)) from observations of the high-rcdshift 
universe may also manifest itself in the properties of the very local Hubble flow. 
We introduce the concept of the critical distance tq where the repulsive force of 
dark energy starts to dominate over the gravity of a mass concentration. For the 
Local Group tq is about 1.5 Mpc. Intriguingly, at the same distance 1.5 Mpc the 
linear and very "cold" Hubble flow emerges, with about the global Hubble constant. 
We also consider the critical epoch tQ, when the DE antigravity began to dominate 
over the local matter gravity for a galaxy which at the present epoch is in the 
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local DE dominated region. Our main result is that the homogeneous dark energy 
component, revealed by SNIa observations, resolves the old confrontation between 
the local Hubble flow and local highly non-uniform, fractal matter distribution. It 
explains why the Hubble law starts on the outskirts of the Local Group, with the 
same Hubble constant as globally and with a remarkably small velocity dispersion. 

1 Introduction 

The Hubble diagram for high redshift Type la Supernovae (Riess et al. 1998, Perlmutter 
et al. 1999) and the Boomerang and MAXIMA- 1 measurements of the first acoustic peak 
location in the angular power spectrum of the cosmic microwave background (CMB) (de 
Bernardis 2000; Jaffe et al. 2000) strongly favor a standard cosmological model having 
the critical density (ft = Q m + Q\ = 1) and a dominant A-like, "dark energy' component 
at the present epoch ps 0.7). In this Letter we argue that the cosmological A 
component is also important for understanding the very local Hubble diagram. Namely, 
there is a recognized riddle in observational cosmology - the linear and quiet Hubble 
flow in the very clumpy local universe (Sandage et al. 1972; Weinberg 1977; Peebles 
1992; Baryshev 1994; Karachentsev & Makarov 1996; Baryshev et al. 1998; Sandage 
1999). 

We show that the Hubble flow emerges there, at a distance of about 1.5 Mpc, where 
the antigravity of the dark energy starts to dominate over the gravity of clumpy matter. 
We suppose that this is not a coincidence but the key to the riddle of the local Hubble 
law. 

2 The enigma of the local Hubble flow 

2.1 Formulation of the problem 

In the standard cosmology the dynamics described by the Hubble law, V = H x r, is a 
strict consequence of a uniform distribution of self-gravitating matter (Robertson 1955; 
Peebles 1993). Hence the predicted linear velocity field is only valid for scales where 
the universe is uniform. However, a very puzzling fact, long ago noted by Sandage, 
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Tammann & Hardy (1972; STH), is that Hubble discovered his law in the distance 
interval 1 — 20 Mpc, where the galaxies are very clumpy distributed. Indeed, this was 
deep inside a cell of uniformity, the size of which is still debated, but there are modern 
estimates from 30 Mpc up to 200 Mpc (Sylos Labini, Montuori & Pietronero 1998; Wu, 
Lahav & Rees 1999) Inside such cells the galaxies are fractally distributed, being clumpy 
at all smaller scales. This means that "in fact, we would not expect any neat relation of 
proportionality between velocity and distance for these [closeby] . . . galaxies" ( Weinberg 
1977). 

On the other hand, recent studies of the local volume (< 10 Mpc), based on accurate 
distances, confirm a "cold" linear Hubble flow, with about the global Hubble constant Hq 
(Sandage 1986; Karachentsev & Makarov 1996; Teerikorpi 1997; Giovanelli et al. 1999: 
Ekholm et al. 1999). Sandage (1999) recently expressed this surprising situation, as the 
"extremely local rate [of expansion] is the same as the global rate to better than 10%", 
and an " explanation of why the local expansion field is so noiseless remains a mystery" . 
Indeed, the velocity dispersion around the local Hubble law is very low (<jy < 50 — 70 
km/s) (Sandage 1986; Karachentsev & Makarov 1996; Ekholm et al. 2000; in the last 
work Cepheid distances gave « 40 km/s). That here is a real problem, has been also 
shown with N-body simulations (Governato et al. 1997): the expected velocity scatter 
close to galactic clumps like our LG is 150 — 700 km/s, depending on the CDM model. 

One may thus divide the enigma of the local Hubble flow into the following Problems: 

(1) The linear velocity law appears at distances (1.5 Mpc) which are about 1 percent 
of the scale where the Universe may finally appear uniform. Why does our highly non- 
uniform closeby environment expand as if it were uniform? Why is the rate of expansion 
the same locally and globally, i.e. deep inside and well outside of the cell of uniformity? 

(2) What makes the local Hubble flow so quiet? Not only does the Hubble law exist, 
but it has a remarkably small scatter. 

2.2 Previous attempts to explain the Hubble enigma 

In their seminal paper, STH discussed the hierarchical galaxy distribution suggested by 
de Vaucouleurs (1971). Wertz (1971) and Haggerty & Wertz (1971) had made calcula- 
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tions on how the Hubble law deviates from linearity inside such a structure. STH did 
not find the predicted non-linearity and rejected the hierarchical model. However, they 
emphasized that the co-existence of the linear Hubble law and the large local clouds of 
inhomogeneous matter still is a dilemma. STH proposed two possible solutions. The lin- 
early expanding universe may be inhomogeneous if its mean density is very low, << p cr it- 
Or it may be filled by a totally uniform dark matter. 

Modern redshift surveys, together with accurate measurements of the local velocity 
field, have sharpened this "Hubble - de Vaucouleurs" paradox (Baryshev et al. 1998). 
Indeed, the 3-dimensional galaxy maps at scales from 0.1 Mpc up to 20-200 Mpc are 
well described by a scale invariant fractal distribution (Davis & Peebles 1983; Pietronero 
1987; Peebles 1993; Sylos Labini, Montuori & Pietronero 1998; Teerikorpi et al. 1998; 
Wu, Lahav & Rees 1999). Q Using new data on large scale structure and on the Hubble 
law, the solutions mentioned by STH were studied in Baryshev et al. (1998). ^,From the 
linear theory of the gravitational growth of density fluctuations (without dark energy), 
they concluded that either Q m < 0.01 or for the uniform dark matter Q^ark — 0.99. 

Nevertheless, these explanations are not sufficient, if considered together with recent 
cosmological observations which give £l m ~ 0.3. Also there is no direct evidence for a 
high-density uniform dark matter. Hence the riddle still exists. 

3 Dark energy and dynamical model 

3.1 Vacuum, quintessence, dark energy 

Recent detection of the dark energy component provide a natural candidate for a high- 
density truly uniform background, giving a novel possibility for resolving the Hubble 
enigma. 

In the cosmological theory there are several kinds of dark energy with a positive 

energy density {cq = pqc 2 > 0) and negative pressure (jjq < 0). The dark energy 

or quintessence is a common name for: the classical Einstein's cosmological constant 
1 Note that Mandelbrot (1982) replaced the old concept of hierarchical galaxy distribution by the 
more adequate concept of stochastic fractals. 
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A, time- variable cosmological constant A(t), cosmological vacuum, effective scalar field 
and other forms of exotic matter with negative pressure (see reviews by Bahcall et al. 
1999; Sahni & Starobinsky 2000). Dark energy is now much discussed and used for 
interpretation of observational cosmological tests (Wang et al. 2000; Podariu & Ratra 
2000). The main reason to consider evolving DE is the "cosmic coincidence": Why is 
the rapidly decreasing energy density of matter at the present epoch so close to that of 
the unchanging vacuum density? 

The most striking property of the dark energy having the equation of state pq = weq , 
with w € [—1,0), is that its gravitating mass 

Mq = ^-(l + 3w) PQ r 3 (1) 

is negative for w < —1/3. This produces cosmological " antigravity" or a repulsion 
which accelerates the expansion of the universe. E.g. for Einstein's A = 8-itGpa/c 2 
(cosmological vacuum with w = — 1), the gravitating mass is Ma = — 87rpAr 3 /3. 

3.2 Dynamical model and the critical distance r Q 

We consider the dark energy component uniformly spreading everywhere, and existing 
alongside with luminous and dark matter in the immediate environment of our Galaxy. 
Hence, the local dynamics is determined by the competition of the gravity of dark 
(and luminous) matter and the antigravity of the DE. The mass-to-luminosity ratios in 
systems of different scales suggest that M/L remains constant for r > 0.5 Mpc (Bahcall, 
Lubin, & Dorman 1995), hence dark matter should be distributed like luminous matter 
on such scales. The correlation analysis of the space distribution of galaxies shows 
that the mean density of luminous matter (and hence, of dark matter) decreases with 
increasing scale so that the matter mass increases as M m (r) oc r D , where D, the fractal 
dimension, is between 1 and 2 at least on scales up to 20 Mpc (Wu, Lahav & Rees 1999 
). As its density is constant (Mg(r) oc r 3 ), the DE component starts to dominate after 
some distance tq. 

For a simple estimate of the scale beyond which vacuum or quintessence dominates, 
we consider two models that may imitate the environment of the Local Group, where the 



5 



local Hubble law is accurately observed. The first model is a generalization of Sandage's 
(1986) point-mass model, where a mass M is placed on the dark energy background with 
density pQ. In the second model, a spherical distribution of mass M m (r) is placed on 
the DE background. 

The dynamics of expansion is described by the equation of motion which is the (1,1)- 
component of Einstein's field equations under the assumptions of spherical symmetry, 
dust-like matter, and cosmological dark energy given by the energy-momentum (EM) 
tensor Tj: = pqc 2 diag(l, —w, —w, —w): 

r = -GM eff /r 2 ; M eff = M m (r) + M Q {r) . (2) 

Here M m (r) = M, or M m (r) = M*(r/r*) D , respectively for the first and second model, 
M Q is given by Eq.@. 

For the point-mass the critical distance tq corresponding to r = 0, i.e. the distance 
where DE antigravity compensates matter gravity, is defined by: 

r>r Q = (3M/(4vru}p Q )) 1 / 3 (3) 

where w = — 3w — 1. For w > —1/3 there is no repulsive force and tq does not exist. 
For the second model the dark energy term dominates dynamically at distances 

D_ _j_ 3 _J_ 

r>r Q = (r,) s-u (M^s-d (— — )s-d (4) 

3.3 The temporal behavior of tq and critical time tQ 

In the case of cosmological constant a galaxy which now is on the border of the vacuum 
dominated region around a mass clump, previously lied inside the gravity dominated 
sphere. In terms of comoving coordinates, tq shifts outwards for increasing redshift. 
This prompts us to ask how tq behaves if the DE density changes, as in quintessence 
models, where rq may still shift outwards in comoving coordinates, but slower than for 
the constant vacuum. 

We illustrate the shift of tq in time for three types of DE models: 1) Einstein's cos- 
mological constant A, 2) the quintessence with w = —2/3, and 3) the case of "coherent" 
evolution of a quintessence (pq = kphm, where p\ im is the uniform matter component) 
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during the late history of the universe, with w = —2/3, k = 1. The condition that the 
covariant divergence of the total EM tensor (quintessence + matter) is zero implies that 
p = — 3(/9+ ^)d/a, where p = pq + p m and p = Pq for dust-like matter (p m = 0). Hence 
the DE density behaves as 

n( l-\-k-\-wk \ 

pqoca ^ i+* > (5) 

For 1/k = (no homogeneous matter) Eq.([|) gives the usual (non-coherent) quintessence 
behavior pq oc a~ s ( 1+w \ 

Let us consider the ratio between two metric distances r ga [/rQ, where r ga \ is the 
distance to a galaxy which now is in the DE dominated region and takes part in the 
Hubble flow. The distance r ga i is simply proportional to a. The behavior of the critical 
distance rn depends on the DE and matter models. E.g. for the point-mass model the 
ratio is r ga i/rQ = a( _m /( 1 + 1 / fc )) ) which for 1/k = gives a~ w , corresponding either to 
Einstein's A (w = —1) or to the quintessence (for w < —1/3). 

We define a characteristic time for the antigravity dominance. It is counted from 
the critical epoch ig, corresponding to the critical scale factor ag = a(t = tq) when DE 
antigravity starts to dominate over the gravity of matter concentration component for 
a galaxy which presently [t = to) is a t the distance r = r ga i = 2tq, i.e. well within the 
DE dominated region. E.g. for Sandage's point-mass model thus defined critical scale 
factor is 

a Q = 2 ( 1 + 1 A)/- (6) 

With 1/k = this gives ag = 0.5 for w = —1 and ag = 0.354 for w = —2/3, while with 
k = l,w = —2/3 the critical epoch is at ag = 0.125. This illustrates the difference in 
critical times for different quintessence models. 



4 Specific properties of dark energy dominated dynamics 
4.1 Structure evolution in DE dominated regions 

The theory of structure formation in the general case of an evolving DE component 
is still under construction, though there are some results for particular cases (Sahni &; 
Starobinsky 2000; Fabris &; Goncalves 2000). If during the structure formation period 
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most of the energy of the universe produces antigravity and hence resists gravitational 
collapse, it is impossible for the structure to grow at all. On the other hand, if antigravity 
dominates only a short period, as in the case of the cosmological constant, then even 
for ~ 1 it has only a slight effect on the infall velocities around growing structures 
(Peebles 1984; Lahav et al. 1991). This is clear from Eq.(|2|) which gives for the ratio 
of DE and point-mass accelerations gQ/g m ct a~ 3w , hence for w = —1 the dynamical 
influence of the constant A-component decreases very fast, while for w = —2/3 the 
antigravity acts longer. 

Another important consequence of the DE component is that within A dominated 
regions initial peculiar velocities decrease due to adiabatic cooling. The linear analysis of 
gravitational instability, following Zeldovich (1965), shows that only decreasing modes 
exist in the vacuum-dominated region (Chernin, Teerikorpi & Baryshev 2000). This 
property of the accelerated universe is important for understanding the local Hubble 
flow. Our simple dynamical models led to the notion of the distance tq to the border 
between the matter and DE dominated regions. If a galaxy now beyond r = tq also in 
the past was long enough in the DE dominated region, this would readily explain the 
smooth Hubble law at such distances. Then initial peculiar velocities induced by the 
mass of the Local Group (and similar more distant groups) have decreased because of 
the adiabatic cooling (SV(now) = 5V(z)a(z)). 

Therefore a compromise situation is possible when the time- variable DE is dominant 
during a period, sufficient for the structure growth in early times and for the cooling of 
velocity dispersion in late times. This means that such DE dominated regions between 
mass concentrations are "pacific oceans" where the Hubble law appears and the global 
Hubble constant may be measured even locally. 

4.2 Restoration of the Hubble flow in volumes with bulk motion 

We referred to two properties of the cosmological vacuum, namely uniformity and anti- 
gravity. Let us now consider an important third property of the vacuum. It is clear that 
the medium with the equation of state p = —pcP described by the cosmological constant 
has EM tensor which is proportional to the metric tensor = agf. Hence it has the 
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remarkable mechanical property of vacuum: rest and motion can not be discriminated 
relative to it. In other words, any matter motion is co-moving to this medium. 

When the vacuum dominates over the self-gravity of matter, matter masses like 
galaxies move as test particles in the antigravity of the vacuum. The vacuum accelerates 
the motions tending to form a regular kinematic pattern with a linear velocity field. The 
expansion rate in such a flow depends only on the value of the vacuum density. 

These considerations can be applied to a local vacuum-dominated area, no matter 
how fast is the bulk motion of the area against the CMB. They are valid as well for 
the large-scale matter distribution; in both cases, the flow will be linear and the rate of 
expansion will be the same on any scale. And the flow will also be stable against velocity 
or density perturbations. 

5 The critical distance vq for the Local Group 

Let us calculate the critical distance for the Local Group in the case of the point-mass 
model and the cosmological vacuum (w = — 1 and pQ = p\ is constant in time). 

The most recent estimate for the mass of the LG gives Mlg = 2 x 1O 12 M (van 
den Bergh 1999). The vacuum density from SNIa observations is p\ ~ 0.7p cr u = 4.7 x 
lCT 30 /ig g/cm 3 . Then the distance r\, where the vacuum starts to dominate, is ~ 1.5 
Mpc in the point-mass model. If one changes Mlg or Pa by a factor of two, the critical 
distance is changed by 26 percent. If one chooses the mass distribution model with 
fractal dimension D = 1 and M* = 2 x 10 12 M Q at r* = 1 Mpc, then Eq. (4) gives 
r A = 1.8 Mpc. 

Hence, r\ is robustly put into the range from 1 to 2 Mpc. This value is surprisingly 
close to the distance where the Hubble law emerges, which is 1.5 Mpc (Sandage 1986). 
Is this just a coincidence? Not at all, we suggest that this is rather a major feature of 
the local matter flow. Indeed, the dynamical dominance of the dark energy at 1-2 Mpc 
and beyond means that the homogeneous vacuum provides the dynamical background 
for the Hubble flow of matter, on these distances. 
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6 The local and global values of the Hubble constant 



Now we can understand Problem I - the near-equality of the local and global Hubble 
constants - if the undistorted Hubble law is the signature of the dynamical dominance 
of the uniform DE component. In regions, which have for a sufficient time been DE- 
dominated, the rate of expansion - in the first (and main!) approximation - is expressed 
via the DE density: Hq rs (87tG/?q/3) 1//2 . The lumpy, fractally distributed matter gives 
only a small correction to this formula for the global Hubble constant, because its average 
density is much less than the DE density. The very interesting aspect of the dominating 
uniform DE component is that it allows us to predict the value of the Hubble constant 
from a local estimate of the dark energy density. 

If we identify the distance r s t ar t where the Hubble law starts, with the critical radius 
tq for our Local Group, then the predicted value of Hq is 

H(predicted) » ^^)^ (r start )-^ (7) 
to 

where M^q is the mass of the Local Group. With the "standard" values M = 2 x 1O 12 M 
and r s tart = 1-5 Mpc, the Eq.(^) gives for Hq values from 53 km/s/Mpc (w = —1) to 
74 km/s/Mpc (w = —2/3), covering the values obtained for Hq by different research 
groups. More conservatively, these are upper limits, if in fact tq < r s t ar t- We regard 
this agreement with direct measurements of Hq as evidence for the cosmological vacuum 
and quintessence models. Note that this reasoning for the DE component is based 
on locally measured quantities (^Adj^Q and v S f ar f ) and is independent of the supernovae 
determination of $7 a (which did not depend on the value of Hq). 



7 The quietness of the local Hubble flow 

The possibility to solve Problem II may be illustrated by the following simplified rea- 
soning. At any epoch a(t) and on the scale r one may regard 5V ~ Hr as an upper limit 
to the scatter around the Hubble law. For models considered in Sec. 3. 3 the curvature 
is zero and hence from Eq.([|) it follows that the Hubble constant H(a) depends on the 
scale factor as a~( 3 / 2 ^, where j3 = (1 + k + wk)/(l + k). For example in the case of 
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vacuum, i.e. w = — 1 and 1/k = we have H = const, and in the case of dust matter 
k = o hence H oc a~ 3 / 2 . 

So the scatter SV(a) < HqTq a 1 ~( 3 / 2 ) /3 , where tq is the considered scale at the present 
epoch. E.g. for r = 5 Mpc and H = 60 km/s/Mpc, 6V(a) < 300 a 1 ^ 3 / 2 ^ km/s. This 
is also valid at the critical epoch oq after which the peculiar velocities within a DE 
dominated region start to cool down adiabatically, so that 5V(now) = 5V((1q) ag. 
Hence one expects now 5V(now) < 300 ag - ^ 2 ^ km/s. 

For simple numerical examples we use the critical epoch ciq calculated in Sec. 3. 3. In 
the case of the cosmological constant we have oq = 0.500, hence 8V(now) < 75 km/s. 
For the quintessence with w = —2/3, oq = 0.354 and 5V(now) < 63 km/s, and for the 
case of coherent evolution w = —2/3, k = 1, the critical scale factor clq = 0.125 and now 
one expects to observe SV (now) < 38 km/s. 

Of course, these numbers just illuminate the trend for different DE models. A com- 
plete treatment of the dispersion problem needs high-resolution large- volume simulations 
of structure formation on the evolving DE background. 

8 Conclusions 

We conclude that a cosmological dark energy component allows one to understand the 
emergence of the linear and quiet Hubble flow at the border of the Local Group. The 
coincidence of the tq with the starting distance of the Hubble law is seen to have the 
same roots as the "cosmic coincidence" of A and matter densities. A summary of our 
results is: 

* We introduce the concept of the critical distance tq where the repulsive force starts 
to dominate over the gravity of matter. For the Local Group tq « 1.5 Mpc in the case 
of pQ = 0.7p cr i t . This is close to the distance where the local Hubble flow emerges. We 
make this new cosmic coincidence as a starting point for understanding the cold local 
Hubble flow as reflecting the dark energy dominated region around us. 

* In dark energy dominated regions of the universe ("pacific oceans") the linear 
Hubble law exists due to the homogeneity of the DE component, with the same Hubble 
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constant at first approximation determined by the dark energy density. This resolves 
Problem I that the Hubble law not only starts at the border of the DE dominated region, 
but it appears with the Ho as expected from the dark energy density. 

* The solution of Problem II follows from the fact that a galaxy spending a sufficient 
time in the DE dominated region looses its initial peculiar velocity which cools down 
adiabatically. The mysterious quietness of the local Hubble flow does not look at all 
dramatic in the context of the dark energy driven expansion. 

* Our results support a dark energy component which presently dominates the dy- 
namics of the Universe on scales from the Local Group to the Hubble radius. Even 
the closeby galaxy universe is now seen as a cosmic laboratory where all the physical 
ingredients of the universe: luminous matter, dark matter, and the dark energy, may be 
detected. Hence, the very local volume is extremely important for the study of global 
properties of the universe. Cosmology starts immediately beyond the border of the Local 
Group! 

In other environments the DE dominance may start from smaller or larger distances. 
So, around the Coma cluster, the vacuum starts dominating from about 20 Mpc. For a 
cosmologist in Coma, cosmology begins around such a distance. In such DE dominated 
regions the quiet Hubble flow is produced by the same vacuum density and hence the 
Hubble constant is the same. 

One now realizes that when Hubble found his linear redshift law in the closeby 
lumpy environment, he actually saw the influence of the mysterious uniform dark energy. 
Though the value of the vacuum or DE density is still poorly understood in cosmological 
theory, the observations make one agree with the old words by Lemaitre who was the 
first to associate A with the vacuum having positive energy density and negative pressure 
leading to a cosmological repulsion: "Everything happens as though the energy in vacuo 
were different from zero." (Seitter & Duemmler 1989). 
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